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We quantify the scaling of magnetic energy density in the inertial range of solar wind turbulence 
seen in-situ at 1AU with respect to solar activity. At solar maximum, when the coronal magnetic 
field is dynamic and topologically complex, we find self-similar scaling in the solar wind, whereas at 
solar mimimum, when the coronal fields are more ordered, we find multifractality. This quantifies 
the solar wind signature that is of direct coronal origin, and distinguishes it from that of local MHD 
turbulence, with quantitative implications for coronal heating of the solar wind. 
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The interplanetary solar wind exhibits fluctuations 
characteristic of Magnetohydrodynamic (MHD) turbu- 
lence evolving in the presence of structures of coronal ori- 
gin. In-situ spacecraft observations of plasma parameters 
are at minute (or below) resolution for intervals spanning 
the solar cycle, and provide a large number of samples 
for statistical studies. These reveal a magnetic Reynolds 
number ~ 10 5 [l[ and power spectra with a clear inertial 
range over several orders of magnitude characterised by a 
power law Kolmogorov exponent of ~ —5/3. Quantifying 
the properties of fluctuations in the solar wind can thus 
provide insights into MHD turbulence and also inform 
our understanding of coronal processes and ultimately, 
the mechanisms by which the solar wind is heated. Quan- 
tifying these fluctuations is also central to understanding 
the propagation of cosmic rays in the heliosphere Q • 

Coronal heating mechanisms are studied in terms of 
the scaling properties of coronal structures 0, Q , heat- 
ing rates [5J1 and diffusion via random walks of magnetic 
field lines [6[, all of which suggest self-similar processes. 
The solar wind is also studied in-situ to infer informa- 
tion pertaining to coronal processes. Large scale coherent 
structures of solar origin, such as CMEs, can be directly 
identified in the solar wind. At frequencies below the 
'Kolmogorov- like' inertial range, the solar wind exhibits 
an energy containing range which shows ~ 1// scaling 
0- Solar flares show scale invariance in their en- 
ergy release statistics over several orders of magnitude 
Q which has been discussed in terms of Self-Organized 
Criticality (SOC) 0,[lH. Within the inertial range, the 
observed solar wind magnetic fluctuations are principally 
Alfvcnic in character with asymmetric propagation anti- 
sunwards [121 ]. In-situ plasma parameters which directly 
relate to cascade theories of ideal incompressible MHD 
turbulence, such as velocity, magnetic field, and the El- 
sasser variables have thus been extensively studied in the 
solar wind ([l3| and refs. therein). These show multi- 
fractal scaling in their higher order moments consistent 



with intermittent turbulence 1J, [15| . Intriguingly, the 
magnetic energy density B 2 and number density p show 
approximately self-similar scaling in the inertial range 
lg,liZ|- These parameters are insensitive to Alfvenicity, 
and do not relate directly to MHD cascade theories. 

In this Letter we quantify the scaling seen in B 2 in the 
inertial range of solar wind turbulence with respect to 
coronal structure and dynamics. We employ a recently 
developed technique [18| that sensitively distinguishes 
between self-similarity and multifractality in timeseries. 
This will allow us to distinguish and quantify the solar 
wind signature that is of direct coronal origin from that of 
local MHD turbulence, with quantitative implications for 
our understanding of coronal heating of the solar wind. 

The WIND and ACE spacecraft spend extended inter- 
vals at ~1 AU in the ecliptic and provide in-situ magnetic 
field observations of the solar wind over extended peri- 
ods covering all phases of the solar cycle. We focus on 
a comparison between solar maximum when the coronal 
structure is highly variable with topologically complex 
magnetic structure, with that at solar minimum when the 
coronal magnetic structure is highly ordered. The most 
magnetically ordered region of the corona is at the poles 
at solar minimum and observations of the correspond- 
ing quiet, fast solar wind are provided by the ULYSSES 
spacecraft. The four data sets, [23], that we consider here 
are then a.) WIND 60 seconds averaged MFI data at the 
solar maximum year of 2000 and 6.) at the solar min- 
imum year of 1996; c.) ACE 64 seconds averaged MFI 
data for the year 2000; and d.) ULYSSES 60 seconds 
averaged VHM/FGM data for July and August 1995. 
Data sets a-c consist of ~ 4.5 x 10 5 points; and d con- 
sists of ~ 8.5 x 10 4 data points. Intervals corresponding 
to magnetospheric bow shock crossings for WIND were 
removed by comparison with [28|. The ACE spacecraft 
orbits around the Earth-Sun LI point and the ULYSSES 
data was obtained for the North polar pass of 1995. All 
of the above intervals show a ~ —5/3 power law scaling 
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inertial range in the power spectra of \B\ over several 
decades which is indicative of a well developed turbulent 
fluid. 

We can access the statistical scaling properties of a 
timeseries by constructing differences y(t,r) — \B(t + 
t)| 2 — \B(t)\ 2 on all available time intervals r. The statis- 
tical scaling with r can be seen in the structure functions 
of order m which follows that of the moments of the PDF 
of y, P(y,r): 

/CO 
\y\ m P(y,r)dy , (1) 
-oo 

where (} indicate ensemble averaging over t. Statistical 
self-similarity implies that any PDF at scale r can be 
collapsed onto a unique single variable PDF V s : 

P(y,T)=T- H V s (T- H y) , (2) 

where H is the Hurst exponent. Equation @ implies 
that the increments y are self-affine i.e. they obey the 
statistical scaling equality yibr) = b H y(r) , such that the 
structure functions will scale with r as 

S m ( T ) =t^S?(1) , (3) 

For the special case of a statistically self-similar (ran- 
dom fractal) process, £( TO ) — Hm. This scaling with 
H = 1/3 is characteristic of Kolmogorov's 1941 theory 
of turbulence and intermittency corrections to this 
are modelled by quadratic and concave C( m ) (multifrac- 
tals) [2{j. A difficulty that can arise in the experimen- 
tal determination of the £( m ) is that for a finite length 
timeseries, the integral ([1) is not sampled over the range 
(— oo, +oo); the outlying measured values of y determine 
the limits. In the case of a heavy-tailed PDF the higher 
order moments (larger m) can yield a ((m) that deviates 
strongly from the scaling of P(y, r) in ((2|) [l8| (hereafter 
KCH) . An operational solution to this was demonstrated 
in KCH for a self-similar process. Essentially one sys- 
tematically excludes a minimal percentage of the outly- 
ing events y from the integral in {T]) so that the statistics 
of the PDF tails become well sampled and the integral 
(TTJ) yields a r dependence with the correct scaling of the 
self-similar process This method is sensitive in dis- 
tinguishing self-affine scaling from weak multifractality. 
We illustrate this with two reference models: the first 
of which is manifestly self-similar, an a-stable Levy pro- 
cess of index a = 1.0 (H = 1/a) (HI; and the second, 
manifestly multifractal, i.e. a p-model (2ll | with p = 0.6. 
These synthetic data sets each consist of 10 6 data points. 
Figure [T] shows plots of the exponents C( TO ) Vs. m ob- 
tained from (J3]) by computing the gradients of log S m (r) 
for (a) the Levy process and (b) the multifractal model 
respectively. The exponents have been recomputed 
as outlying data points are successively removed, and we 
can see that removing a small fraction, ~ 0.001% of the 



data leads to a large change in the computed C(p). A 
reliable estimate of the exponents from the data requires 
rapid convergence to robust values; shown in KCH to be 
a property of self-affine timeseries. We can see this be- 
haviour in the Levy model which quickly converges to 
linear dependence of with p as expected. 
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Figure 1: C( m ) Vs. m plots for a.) a — 1.0 symmetric Levy 
process and b.) p = 0.6 p-model process. 

The multifractal p-model only begins to approach lin- 
earity after ~ 3% of the data is excluded. This apparent 
linearity in the p-model is actually a divergence in the 
values of the C(p)- We see this behaviour if we plot the 
value of one of the exponents from Figure [T] versus the 
percentage of points removed. This is shown for £(2) for 
the Levy process (upper panel) and the p-model (lower 
panel) in Figure [5] As we successively exclude outlying 
data points, the self-affine Levy process quickly reaches 
a constant value for ((2) — 2/a — 2.0 whereas for the 
multifractal, the £(2) exponent shows a continuing sec- 
ular drift. Importantly, successively removing outlying 
data points does not convert the multifractal p-model 
timeseries into a self-affine process. In addition, a plot 
of versus p (Figure [lj is not sufficient to distinguish 
self-affine from multifractal behaviour, one also needs to 
examine the convergence properties of the exponents as 
outlying points are successively removed, as shown in 
Figure [U 

We now turn to the analysis of solar wind data. In 
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Figure 2: Exponent of the second moment ("(2) Vs. the per- 
centage of points excluded for a.) the Levy model and b.) 
p-model. 
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Figure 3: Exponent of the second moment £(2) Vs. the per- 
centage of points excluded for a.) WIND and ACE at solar 
maximum and b.) WIND at solar minimum. 
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Figure 4: Exponent of the second moment £(2) Vs. the per- 
centage of points excluded for ULYSSES at solar minimum. 



Figures [3^l and b we plot £(2) versus the percentage of 
points removed in B 2 for intervals at solar maximum and 
minimum respectively. The C values for these plots were 
obtained from an identified scaling range which spanned 
from ~ 5.2 minutes to ~ 2.7 hours (see e.g. [lfjl . |22| 
). Comparison of these plots with Figure [2] strongly 
suggests that at solar maximum, the magnetic energy 
density is self-affine; we can clearly identify a plateau 
with a H = C(2)/2 value of H ~ 0.44 ± 0.02 for WIND 
and H ~ 0.45 ± 0.01 for ACE. At solar minimum, there 
is no clear plateau and the behaviour is reminiscent of 
the multifractal p-model. We have thus differentiated 
the distinct scaling behaviour at solar maximum and so- 
lar minimum. Intriguingly, it is at solar maximum that 
we see self-similar behaviour; whereas at solar minimum 
the timeseries resembles a multifractal, reminiscent of in- 
termittent turbulence. Since the corona is complex and 
highly structured at solar maximum, this is highly sug- 
gestive that this self-similar signature in B 2 is related 
to coronal structure and dynamics rather than to local 
turbulence. 

We can test this conjecture by considering observations 
of the solar wind where the coronal structure is maxi- 
mally ordered. We repeat the above analysis on a two 
month interval of ULYSSES data during solar minimum. 
The resulting plot of ((2) versus percentage of points 
excluded is shown in Figure [4] This plot again does not 
support self-affine scaling and is reminiscent of that of the 
p-model, strengthening previous results 15,23|. Clearly, 
the behaviour of B 2 in the solar wind originating from 
a corona dominated by ordered open field lines is not 
self-affine. The appearance of fractal versus multifractal 
behaviour in B 2 is not a strong discriminator of vari- 
ability in the average solar wind speed per se. We see 
multifractal scaling both in the ecliptic at minimum in 
an interval that contains periods of alternating high and 
low speed streams, and at the poles, where the average 
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speed is fast and uniform. We see fractal scaling in the 
ecliptic at maximum, where the average speed also alter- 
nates between fast and slow streams. Previous studies 
(23| have shown a variation with latitude and solar cycle 
of the level of multifractality of components of magnetic 
field. This may be related to the signature of the level of 
complexity in coronal magnetic structure which we have 
identified in B 2 within the inertial range of turbulence, 
but may also simply reflect a correlation with average so- 
lar wind speed. We have also verified that \B\ does not 
show evidence of self-similarity for the intervals chosen 
for our study. More specifically \B\ exhibits multifractal 
behaviour. This confirms the earlier results of Hnat et. 

al. m. 

The corona contains many long-lived structures which 
extend far out into the solar system mediated by the 
interplanetary solar wind At solar maximum these 
structures show a high degree of topological complexity. 
One model for these structures and their propagation is 
as a random walk or braiding of magnetic field lines with 
a measurable diffusion coefficient [3, IS ___} • A diffusion 
process such as this intrinsically generates self-similar 
scaling, and may in a straightforward manner account for 
that shown here in B 2 at solar maximum. Alternatively, 
the relevant process may be that of reconnection in the 
complex magnetic structure of the emerging coronal flux. 



high-frequency self-similarity reported here may suggest 
further processes responsible for coronal heating. 
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Models for this include SOC based random networks [11 1 
which again imply self-similar scaling. Our quantitative 
determination of the Hurst exponent H ~ 0.45 of the 
sclf-afime scaling seen in the solar wind provides a strong 
constraint to these models. 

The PDF resulting from such a self-similar process can 
be captured by a solution to a generalized Fokker-Planck 
equation (FPE) with power law scaling of the transport 
coefficients Intriguingly, the associated Langevin 

equation transforms nonlinearly into that for a constant 
diffusion coefficient. The transformation may be equiva- 
lent to introducing a diffusion process with constant dif- 
fusion coefficient, on a space with non-Euclidean, self- 
similar, fractal geometry This may provide a quanti- 
tative basis for models of transport of initially random 
fractal fields (the coronal source) in a turbulent flow (the 
solar wind). At solar minimum we see quite a differ- 
ent picture. Here the corona is topologically well or- 
dered magnetically. Thus in this case the fluctuations 
in B 2 are dominated by the evolving turbulence of the 
interplanetary solar wind which is well known to exhibit 
multifractal behaviour. Intriguingly, this self-affine sig- 
nature quantified here in B 2 extends over the ~ —5/3 
exponent inertial range seen in the solar wind. This is at 
higher frequencies than the ~ 1// behaviour previously 
identified as a coronal signature in the solar wind [8J. 
Although models involving reconnection and flares and 
nanoflares have been proposed estimates of the to- 
tal energy contained in such structures falls significantly 
short of that required for coronal heating [26j. Thus the 
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